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O ■ Abstract 

o : 

'""'■ We present Nanolmprint lithography experiments and modeling of thin films of block copolymers 

^ . (BCP). The Nanolmprint technique is found to be an efficient tool not only to align lamellar phases 

m ; 

ff^ , perpendicularly to the substrate, but also to get rid of in-plane defects over distances much larger 

!>■ 

Tlj" ' than the natural lamellar periodicity. The modeling relies on self-consistent field calculations done 

t^ ■ 

f^ ' in two- and three-dimensions, and is found to be in good agreement with the experiments. It 

o: 

T— I , also offers some insight on the Nanolmprint lithography setup and on the conditions required to 

. !^ ' perfectly ordered BCP lamellae. 
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I. INTRODUCTION 

One of the main challenges of contemporary design of microchips is to find affordable 
techniques of patterning silicon wafers at the nanoscopic level [1] . To address this question, 
self- assembling block copolymers (BCP) have been suggested as potent candidates to provide 
patterns for nanolithography [2, 3]. The aim is to produce thin and structured BCP films 
with patterns such as lamellae, cylinders and spheres [4] that can then be transferred to a 
substrate. 

When thin films of BCP are cast on a surface, they self-assemble into one of several possi- 
ble nanostructures having a specific orientation with respect to the substrate. In particular, 
by adjusting the surface interactions and film thickness, it is possible to produce lamellar 
and cylindrical phases in an orientation perpendicular [5, 6] to the substrate. However, while 
lamellar and cylindrical phases can be perpendicularly aligned on a large scale [7, 8], the 
ever-remaining challenge for micro-electronic applications is to find affordable and efficient 
techniques for in-plane organization, with minimal amount of defects. This will allow pro- 
ducing devices that are hundreds of micrometers in size, where precise spatial accessibility 
is required. 

Several attempts have been made to address this challenge. They include, among others, 
chemical patterning of the substrate by e-beam lithography [9, 10] or by prepositioning of 
another copolymer layer [11], and graphoepitaxy [12], where an artificial surface topogra- 
phy of grooves separated by walls is created on the substrate. Due to such topographical 
constraints, ordered regions of BCP are obtained over length scales of micrometers [12-15]. 

A more recent technique addressing the same issue is the nanoimprint lithography 
(NIL) [16-19], and it has potential advantages in terms of cost and simplicity. It uses sur- 
face micrometer-sized structural features of a reusable mold made by standard lithographic 
techniques to guide the self-assembly of the BCP at the nanometer scale. In Ref. [18] a 
cylindrical phase oriented perpendicularly to the substrate was imprinted by NIL where the 
NIL patterns are about twice or three times larger than the BCP period. The resulting BCP 
phase is aligned perpendicularly, but the imprint procedure is found to induce unwanted de- 
fects. In another attempt [19], alignment of lamellae was produced by NIL, but the lamellae 
were not oriented perpendicular to the substrate. 

Motivated by previous studies, we attempt in this paper to improve on the NIL procedure. 
In particular, we focus on NIL setups that produce perpendicular aligned lamellae with no 
defects, even when the NIL structural features are much larger than the BCP period. 



II. MATERIALS AND METHODS 

The symmetric di-block copolymer polystyrene-b-polymethylmethacrylate (PS52K-b- 
PMMA52K, PDI:1.09) was purchased from Polymer Source Inc and exhibits a lamellar phase 
of period £0= 49 nm in the bulk [10]. PS and PMMA blocks share very similar values 
of surface tensions, in the range of 29.7 mN/m - 29.9 mN/m and 29.9 mN/m - 31 mN/m, 
respectively, at the temperatures used for film annealing (~170°C). The glass transition 
temperature of PS and PMMA is 100°C and 105°C, respectively. 

UV/ozone apparatus was used both for wafer cleaning and octadecyl-trichlorosilane 
(OTS) self-assembly monolayer (SAM) oxidation. The wafers were cleaned by irradiation 
before the silanization treatment. The OTS SAM oxidation was performed as described in 
Ref. [6] in order to reach a precise surface energy. To get a better contrast during AFM 
imaging, selective phase etching of PMMA can be applied just before the imaging. To avoid 
ozone formation that would degrade both the PS and PMMA blocks, the reaction chamber 
was flushed for 30min with nitrogen before and during UV-irradiation. Samples were then 
irradiated for 30 min and subsequently treated with glacial acetic acid during Ih to ensure 
complete removal of degraded PMMA fragments. 

Solutions of 1 wt% of BCP in toluene were prepared and spin coated at 1800 rpm onto 
silanized silicon wafers to produce BCP fllms with thickness of about 40 nm (below i^). 
Subsequently, the samples were annealed in a vacuum oven at 170°C at a pressure of less than 
3 KPa for 1 day. The resulting lamellar phase was examined and found to be perpendicular 
to the neutral substrate but with many in-plane defects [6]. Then, the sample was treated by 
thermal NIL, which consists of embossing the thin BCP fllm, heated above its glass transition 
temperature, by a reusable mold made of a series of grooves as is shown schematically in 
Fig. 1(a). 

Imprint experiments have been carried out at the CEA/LETI clean room in Grenoble, 
France, using a EVG(r)520HE press. The mold is a 4" silicon wafer with topographical 
features made of grooves of tens of nanometers in height with groove width (uh) and inter- 
groove separation (ui) in the range of hundreds of nanometers [Fig. 1(b)]. The mold surface 
is coated with a perfluorinated polymer layer to avoid BCP adhesion. An overall pressure of 
about 0.3 MPa was applied to the BCP fllm, flrst at a temperature of 120°C during 7hours, 
then at 170°C during 60 hours. 

Surface imaging after imposing NIL on the BCP fllm was performed by SEM and AFM. 
SEM imaging was performed with a Hitachi 9300 apparatus, operated at optimal voltages 
of about 500 V. This voltage allows observation of BCP phase organization, without PMMA 
etching, due to a smaller penetration depth of the electron beam and enhanced secondary 
electron emission yield. AFM imaging has been performed using a NanoScope V (Digital 
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FIG. 1: The NIL setup, (a) A grooved top mold is pressed upon a BCP film oriented perpendicularly to a 
bottom substrate. The lamellae follow the direction of the mold grooves, (b) A cut (side view) through 
the modeled NIL setup consisting of a periodically patterned mold (top surface). The periodicity in the 
x-direction \s d = oji + ojh, with ui and ujh being the finger width and groove width, respectively. The 
BCP film fills the gap between the two surfaces and its thickness varies between Li and L^. 

Instruments). The samples were analyzed using the AFM tapping mode in air and with 
silicon cantilevers of 125 /im in length (Ultrasharp, Micro Masch). Their resonance frequency 
ranges between 265 and 400 KHz, whereas their force constant lies in between 20 - 75 N/m 
and the tip radius of curvature is less than 10 nm. The scan rate was chosen in order to 
obtain the best contrast in phase images and the least deviation between height trace and 
retrace scans. 

III. THE SELF-CONSISTENT FIELD THEORY 



We use self-consistent field (SCF) theory to investigate the behavior of a melt of A-B 
di-block copolymer (BCP) film at nano-patterned surfaces. The BCP film has n polymer 
chains, each having a length N = Na + Nb in terms of the Kuhn length a, which is as- 
sumed, for simplicity, to be the same for the A and B monomers. Hence, the A-monomer 
molar fraction / = Na/N is equal to its volume fraction. In addition, hereafter we con- 
centrate on symmetric di-BCP, A^^ = Nb having / = 0.5. The symmetric BCP yields 
thermodynamically stable lamellar phases of periodicity £o, as temperature is lowered below 
the order-disorder temperature (ODT). We rescale all lengths, r — )■ r/£o, by the natural 
periodicity of the BCP, £o — 4.05/?^, where Rg is the chain radius of gyration R^ = Na^/G. 
Similarly, the curvilinear coordinate along the chain contour, s, is rescaled by A^, yielding 
s — )■ s/N. With these conventions, the free energy for such a BCP film confined between 
the two surfaces is 



- \nQc - — dV [uA{r)(f)A{r) + UB{r)(f)B{r)] 

+ ^y'dV77(r)[0A(r) + 0B(r)-l] (1) 

The system has a total rescaled volume Qi. The closed-packing density (monomer per 
unit area) is po = ^^^j and the Flory-Huggins parameter is x- The dimensionless volume 
fractions of the two components are defined as 0a (^) and (f)B{i^), respectively, whereas ujj{r), 
j = A, B, are the auxiliary fields coupled with (/)j{r), and Qc is the single-chain partition 
function in the presence of the ua and ub fields. The third term represents a surface-energy 
preference, where ua and ub are the short-range interaction parameters of the surface with 
the A and B monomers, respectively. Formally, UA{r) and UB{r) are surface fields and 
have non-zero values only on the surface(s). It is worth noting that due to the above 
rescaling, the variables in Eq. (1) should also be rescaled by: x ~^ ^X) ^ji^^) ^ ^'^j(^) and 
Uj{r) -^ Nuj{r). 

Finally, the last term in Eq. (1) includes a Lagrange multiplier ri[r) introduced to ensure 
the incompressibility condition of the BCP melt: 

0A(r) + 0B(r) = 1 for all r G l^i (2) 

By inserting this condition, Eq. (2), in the surface free energy term of Eq. (1), the integrand 
becomes ua4>a + ub4>b = (yUA — UB)(pA + UB- Hence, u{r) = UA^r) —ub{j-) is the only needed 
surface preference field to be employed hereafter. 

Using the saddle-point approximation, we can obtain a set of self-consistent equations 

^A{r) = X<pB{r) - UAir) + r]{r) (3) 

Wb(^) = X'pA^r) - UBir) + r]{r) (4) 

0A(r) = 7- / dsqAir,s)q\{rJ -s) (5) 

Qc Jo 

<pB{r) = -— dsqBir,s)ql{r,l- f -s) (6) 

Qc Jo 

where the incompressibility condition, Eq. (2), is obeyed, / = Na/N and the single-chain 

free energy Qc is: 

Qc = ^Jd?rq\{rJ) (7) 

The two types of propagators qj (r, s) and g] (r, s) (with j = A, B) are solutions of the 
modified diffusion equation 



with the initial condition qA{r, s=0) =qB{r, s=0) = 1, q\{r, s=0) = qB{r,l—f) and 
g]j(r, s=0) = 5'yi(^)/)- This diffusion equation is solved using reffecting boundary condi- 
tions at the two confining surfaces {z=0 and z=L): dq/dr\z=o = and dq/dr\z=L = 0, while 
periodic boundary conditions are used in the perpendicular direction. 

The top surface with parallel and elongated grooves is modeled as an impenetrable bound- 
ary for the polymer chains. In practice, this is done by assigning a large value to the local 
surface field, UA{r) = ub{t) = —20, for all points inside the grooves. In addition, on the 
whole top surface we include a weak attractive field {u = 0.02), which takes into account 
the preferred surface interaction with the A block. 

We use the same numerical method as that in our previous work [25] to solve these self- 
consistent equations. Fields and densities were calculated on a grid with a spatial resolution 
of A ?^ io/lO or A ?^ ^o/20 for the three- and two-dimensional grids, respectively. First, we 
guess an initial set of values for the {uj{r)} auxiliary fields. Then, through the diffusion 
equations we calculate the propagators, qj and {q]}, and from Eq. (5) and Eq. (6) we 
calculate the monomer volume fractions {<f)j}. Next, a new set of values for {ujj{r)} is 
obtained through Eq. (3) and Eq. (4) and this procedure can be iterated until convergence 
is obtained by some conventional criterion. More details are given in Ref. [25]. 

IV. RESULTS 

Figure 2 shows a top view of SEM and AFM images of a PS52K-b-PMMA52K thin film after 
nanoimprinting with a mold made of a series of parallel grooves [Fig. 1]. From the images, it 
can be seen that the lamellar phase is indeed oriented perpendicular to the bottom substrate. 
Moreover, in Fig. 2(a) the SEM image indicates that the thicker sections of the BCP film 
exhibit lamellae that are nicely oriented in a parallel fashion along the groove vertical walls 
(the y-direction), while the thinner sections of the film are too thin to be visualized by 
the SEM technique. The same in-plane ordered lamellae are shown as an enlargement in 
Fig. 2(b), whereas in Fig. 2(c) the region close to the groove vertical wall is visualized by 
an AFM technique. The AFM image indicates a different in-plane ordering of the thick and 
thin sections of the BCP film. While the thick section (right) shows well ordered lamellae of 
periodicity very close to the bulk one, io, which are aligned by the wall, the lamellae of the 
thin section (left) of the film are less oriented and contain some defects. This observation 
can be attributed to the preferred interaction of the mold vertical walls with one of the two 
BCP blocks. The lamellar film to the right of the wall is thicker in the z-direction (height 
Lh as in Fig. 1(b)) and is in contact with the vertical groove wall, while the film section on 
the left has a smaller height (Li), and the wall has no direct infiuence on it. 

In order to further understand these NIL results, we complement the experiments 
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FIG. 2: (a) Top view of a SEM image of the BCP film after nanoimprinting with a NIL mold with 
Ldi = ijjii = 1.5 /^m. The groove height is L^ — Li = 50 nm. The perpendicularly lying BCP lamellae are 
further oriented along the groove long axis, (b) An enlargement of (a) where several lamellae (left) are 
well ordered in-plane, (c) An AFM top view of an enlarged section of the BCP film close to thin-thick 
boundary (the middle vertical line). To the right, the BCP film is thicker (thickness Lh) and is ordered 
by the groove vertical wall, whereas to the left, the BCP is thinner (thickness Li) and less ordered. 



with self-consistent field (SCF) theory calculations performed on symmetric BCP lamel- 
lar phases [20-24]. The NIL setup is modeled by BCP lamellae having a natural periodicity 
£o, and confined in the z-direction between a fiat and neutral bottom surface {u = 0) at z=0 
and a topographical varying surface (the top mold), as is shown schematically in Fig. 1(b). 
The top surface has the form of elongated grooves (in the ^/-direction) of square cross section 
(in the a;-direction). The down-pointing indentations (fingers) have a width of ui separated 
by grooves with a cross-section width ujh- The BCP film thickness measured with respect 
to the z=0 surface varies between Lh inside the grooves and Li in between the grooves 
[Fig. 1(b)]. The mold preferential interaction toward one of the two blocks is modeled by 



an overall non-zejo surface field, u > 0. 



^ 




FIG. 3: Calculated Sd lamellae in a BCP film of size: L^ x Ly x L^ = W£o x lO^o x OMq and 
Nx = 25.0. The system parameters for the NIL are: u = 0.02, d = AIq, oji = ojh = 24, Lh = 0.64 
and Li = O.SIq. (a) Top-view cut at z = 0.34 and (b) side cut at y = 54 of the simulated NIL 
setup. The perfect in-plane order seen in (a) is accompanied by a perfect perpendicular order in (b). 
For comparison, a top-view cut at z = 0.34 is presented in (c) for a similar system bound by two 
neutral and flat surfaces at z = and z = 0.64- It shows several in-plane defects embedded into the 
perpendicularly oriented lamellar phase. 

We performed several three dimensional (3d) SCF calculations to shed light on the film 
in-plane ordering. The 3d system size is L^ x Ly x L^ = 104 >< 104 ^ 0.64, where L^ is 
set to be less than one BCP periodicity in accord with the experiments (Fig. 2). In our 3d 
SCF calculations, a gradual temperature quench is performed. The starting temperature is 
above the order- disorder temperature (ODT), Nxc — 10.5; hence, inside the BCP disordered 
phase. The temperature is then gradually decreased (or, equivalently, the value of Nx is 
gradually increased) until the system reaches Nx = 25, which is well below the ODT. 

With the simulated NIL mold, a perfect perpendicular lamellar structure is found and is 
shown in Fig. 3(a) as a top-view cut ai z = 0.34, and in 3(b) as a side cut aX y = 54- Note 
that in the calculation. Fig. 3(a), the thinner and thicker sections of the film are equally 
ordered, whereas in the experiments. Fig. 2(c), the in-plane ordering is not as good for the 
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FIG. 4: Calculated 3d lamellae in a BCP film of size: L^ x Ly x L^ = 10^ x 10^ x 1.34. All other 
parameters are the same as in Fig. 3(a). (a) Top-view cut at z = 0.3io and (b) side cut at y = 5£o of the 
simulated NIL setup. The perfect in-plane order seen in (a) is accompanied by a perfect perpendicular 
order in (b). 

thinner section of the film (left side of the figure) as compared with the thicker film section 
on the right side. This lack of in-plane ordering in the experiments might be due to lack 
of equilibration or shear effects at the walls, and more detailed investigations are needed to 
further clarify this point. 

As a further check, we compare the NIL setup with a BCP film confined between two 
neutral and flat surfaces. As can be seen in Fig. 3(c), a perpendicular lamellae structure 
arises with no in-plane ordering when the same gradual temperature quench process is 
repeated for the neutral and flat surfaces. It is clear that there are many in-plane defects in 
this case, in contrast with the NIL mold that induces a strong in-plane ordering. 

We show another 3d calculation in Fig. 4, for which all the parameters are the same 
as in Fig. 3 (a) except that we modify the Lh value to be Lh = 1.3£o- We still obtain the 
perfect perpendicular lamellae for this thicker BCP film, which is consistent with the phase 
diagram of Ref [25]. We conclude that the value of Lh is not strictly limited, but can have 
a range of values that will result in a perfect BCP perpendicular lamellar phase. 

Another important question we addressed is the stability of the perpendicular phase for 
wider grooves (larger Uh)- To further explore this issue we complemented the 3d calculations 
by 2d ones where much larger lateral system sizes (up to L^ x Lz = l.Sio x 80£o) can be 
simulated. Although the two-dimensional calculation cannot infer on the degree of in-plane 
ordering and defects, it is of value because it shows the effects of groove width on the stability 
of the perpendicular BCP lamellae. In Fig. 5(a), an ordered perpendicular phase is found 
for the NIL setup with d = 20£o- This figure should be compared with the 3d calculations 
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FIG. 5: (a) Calculated 2d structures of BCP lamellae in the x—z plane for Nx = 20.0. The system 
parameters are u = 0.02, d = 20^o. ^i = ^h = lO^o. L^ = 1.3io and L/ = 0.3£o- (b) Same as (a) but 
with an increased groove periodicity d = 40^0 with uji = uj^ = 20^o- (c) Top view of a SEM image of 
a PS52K-b-PMMA52K filnn after NIL. For groove width of about 20^o the perpendicular lamellar phase 
persists only near the groove edge (white vertical line on the right). In the region far from the groove 
edge, the lamellae probably change their orientation to one that is parallel to the substrate. The bar 
scale is 0.2 fim. 

as shown in Fig. 3(a) for much narrower grooves oi d = 2C.Q. 

The NIL setup has some limitations as can be seen from our calculations and through 
preliminary experimental results. If we increase the groove periodicity (scaling up both 
ujh and uji), the perfect perpendicular lamellae do not persist. Instead, the film breaks up 
into a mixed morphology, combining a perpendicular phase close to the groove wall with 
a parallel phase induced by the horizontal section of the top surface. This can be seen in 
Fig. 5(b), where all parameters are the same as in 5(a) beside a larger groove periodicity, 
d = 40£o- These findings are in accord with the experimental ones shown in Fig. 5(c), where 
SEM images demonstrate the loss of the perpendicular lamellae, for large enough groove 
periodicity, after distances of about 10£o from the groove wall. 

V. DISCUSSION & CONCLUSIONS 

We address in this article the influence of nano-patterning of surfaces on the orientation 
and alignment of lamellar phases of BCP films, using a NIL technique to produce superior 
perpendicular ordering. The main goal of the NIL is to be able to use surface features 
on length scales larger than the BCP natural periodicity io in order to reduce the cost of 
expensive surface preparation treatments. 
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The effect of the NIL setup on BCP in-plane ordering is clearly demonstrated both in the 
experiments and in the modeling. Without the NIL mold, it is possible to obtain perpen- 
dicular lamellae but with many in-plane defects which cannot be annealed away. However, 
with the NIL mold, wetting of the vertical groove wall induces perfect perpendicular ordering 
with minimal amount of defects, over large lateral distances. 

In our model the perfect perpendicular order is induced by a small surface field u = 0.02, 
chosen to agree with PS-PMMA experimental setup, where the relative difference in the 
surface tension between the two blocks and the surface is about 1%. The value of u is 
limited by two opposing trends. On the one hand, u should be large enough so that one of 
the two blocks would wet the groove vertical wall. On the other hand, u should be small 
enough in order not to interfere with the overall perpendicular ordering. 

Another conclusion from the present study is that the NIL groove periodicity cannot be 
much larger than Iq; otherwise the BCP perpendicular order is lost [see Fig. 5(b) and (c)], 
where the perpendicular order is lost for about d ~ 20£o- 

Our SCF calculations are in good agreement with the presented experiments, and provide 
some insight into the conditions needed to obtain perpendicular lamellar phases with minimal 
amount of in-plane defects. We hope that in the future more detailed three-dimensional 
calculations as well as careful investigations of film rheology will shed more light on the 
fundamental behavior as well as applications of BCP films in presence of nano-patterned 
surfaces. 
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